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Abstract 

Data on J/tp production in inelastic proton-proton, proton-nucleus and nucleus- 
nucleus interactions at 158 A-GeV are analyzed and it is shown that the ratio of 
mean multiplicities of J/tp mesons and pions is the same for all these collisions. 
This observation is difficult to understand within current models of J/ ip production 
in nuclear collisions based on the assumption of hard QCD production of charm 
quarks. 
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1 Introduction 



According to the factorization theorem of perturbative QCD |l[ inclusive cross section of 
a hard process should increase proportionally to A in p+A interactions and proportionaly 
to A 2 in A+A collisions. Models describing J /if) production are built on the basis of this 
prediction (for review see 0). They treat creation of a cc pairs as a hard process and they 
further assume that the initial number of J /if) mesons is proportional to the number of 
charm pairs. Therefore in the absence of medium effects the J / if) cross section is expected 
to increase as A(A 2 ) for p+A (A+A) collisions. 

Experimental results on J / if) production in p+A interactions contradict this naive ex- 
pectation showing an increase of the cross section proportional to A ' 9 . This reduction of 
the ^-dependence is usually explained as being predominantly due to final state interac- 
tions of the J /if) meson (or its premeson state) with nucleons [[|. However models based 
on this picture and parameters fitted to the p+A data in general overpredict recent results 
on J /if> production in central Pb+Pb collisions at 158 A-GeV (for review see Q). This 
reduction of J /if) production in the latter data is usually interpreted as due to interactions 
of J /if) (pre) mesons with surrounding high density matter (ultimately the Quark Gluon 
Plasma) §. 

In high energy A+A collisions (from central S+S to central Pb+Pb) the multiplicity 
of pions and strange hadrons increase proportionally to the number of colliding nucleons 
(participant nucleons) |J. These data and their interpretation in terms of a statistical 
QGP model J7| suggested the question whether a similar dependence may be observed 
for charm and consequently for J /if) production. A simple estimation of the centrality 
dependence of the J/ if) to pion ratio in Pb+Pb collisions indicates that this hypothesis 
may be in fact correct JTJ. 

The aim of this paper is to review available experimental results to obtain information 
concerning the A-dependence of J /if) yield. In particular we study the A-dependence 
of the J /if) to pion ratio using results on proton-proton and nucleus-nucleus interactions 
(Sections 2 and 3) and proton-nucleus interactions (Section 4). We summarize also results 
on the yl-dependence of the open charm yield in p+A interactions (Section 5). 
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2 J/ip Multiplicity in p+p Interactions 



In p+p interactions the J / ip cross section was measured at five different collision energies, 
y/s = 6.8 GeV §, 8.7 GeV @, 19.4 GeV 0, 24.3 GeV and 52 GeV [|J. Most of 



the data are measured for xp > and they are not corrected for the branching ratio to 
the measured decay channel. 

The mean multiplicity of J/ip mesons in full momentum space, (J/ip), is obtained in 
the following way. The xp distribution of J/ip is assumed to be symmetric with respect 
to reflection at xp = 0. The most recent values of the branching ratios J/ip — > fi + + fi~ 



(B w = 0.0601) and J/ip — > e + + e~ (B ee = 0.0602) were used |L3] . The cross sections of 
J/ip production were further divided by the cross sections for inelastic p+p interactions 
at the corresponding collision energy. The latter cross sections were calculated according 
to the parametrization of the experimental data given in [I3|. The resulting values of 
(J/ip) are given in Table 1. 

For comparison with the data on nucleus-nucleus collisions at the SPS the mean J/ip 
multiplicity for p+p interactions at 158 GeV (y/s = 17.3 GeV) is needed. In this energy 
range the energy dependence of the integrated cross section for J/ip production can be 
conveniently parametrized by [|i~4T| : 



where a = 12 and o"o are parameters fitted to the data, a J ^ and rrij/^, are J/ip cross 
section and mass, respectively. This parametrization predicts a decrease of the J /if) yield 
by about 25% when going from y/s = 19.4 GeV to y/s = 17.3 GeV. Thus we can estimate 
(J/ip) to be (2.9 ± 0.5) • 10" 6 at y/s = 17.3 GeV using the measured value of (J/ip) = 
(3.8 ± 0.3) ■ 10~ 6 at y/s = 19.4 GeV (see Table I). The result of the above interpolation, 
shown by the open circle in Fig. 1, agrees with the value estimated in Ref. [14| for 



p+p interactions at 150 GeV using data available at this energy |l(J and an additional 
assumption concerning an unpublished ratio of cross sections. 

The mean multiplicity of negatively charged hadrons (more than 90% are 7r~-mesons) 
in nucleon-nucleon (N+N) interactions at 158 GeV is (h~) = 3.01 ±0.06 [|l5|]. This mean 
multiplicity was calculated as (h~) = ((h~) pp + 2(h~) pn + (h~) nn )/A ) where (h~) pp , (h~) pn 
and (h~) nn are mean multiplicities of negatively charged hadrons for p+p, p+n and n+n 
interactions at 158 A-GeV, respectively [[HJ . 

Taking the value of (J/ip) calculated above we obtain (J/ip) /(h^) = (0.96 + 0.17) • 10 -6 
for N+N interactions at y/s = 17.3 GeV. This ratio is further used for the comparison 
with nucleus-nucleus data. 
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3 J/ip Production in Nucleus— Nucleus Collisions 



The production of J /if) in nucleus-nucleus collisions was measured by the NA38 Col- 
laboration for O+Cu, 0+U and S+U interactions at 200 A-GeV @ and by the NA50 



Collaboration for Pb+Pb interactions at 158 A-GeV fL7| , p~8f] . The procedure which allows 



to calculate the (J /if)) / (h~) ratio from the published data is described below using as an 
example the Pb+Pb results. 

The measured J /if) cross section in minimum bias Pb+Pb collisions is: 

B^ai'f = 21.9 ± 0.2 ± 1.6 fib. 

This cross section refers to the NA50 acceptance < y cm < 1 and — 0.5 < cosOcs < 0.5, 
where y cm is the J /if) rapidity calculated in the cm. system and 6cs is the Collins-Soper 
angle |]I9"|| . In order to get an estimate of the total J /if) cross section we assume that 
the J /if) production for y cm > 1 can be neglected and that the distribution in cosOcs is 
uniform IjTOfl. This leads to a correction factor for the acceptance equal to 4. Based on the 



h+p results at 200 GeV |10| one can estimate that neglecting J /if) yield at y cm > 1 may 
lead to an underestimation of the J /if) multiplicity by less than 30%. A similar conclusion 
is reached when the J /if) rapidity distribution in Pb+Pb collisions is assumed to be similar 
to the rapidity distribution of the mesons measured by the NA49 Collaboration |2(| . 
In addition, the cross section presented by NA50 is corrected here for the branching ratio 
J /if) — > fi + + ji~ (5 W = 0.0601) The cross section for J /if) production resulting from 



the above procedure is: 

^J/V = 1.46 + 0.12 mb, 

where systematic uncertainty of our extrapolation procedure is not included in the quoted 
error. The J / if) multiplicity can be calculated as 

(J/ifj) = = (2.07 ± 0.17) ■ 10~ 4 , 

cr 

where a is the total cross section of inelastic Pb+Pb collisions calculated to be 7040 mb 



using a parametrization of the measured data given in Ref. |2T . 

The results of the NA35 and NA49 Collaborations indicate that the ratio of (h~) to 
the mean number of participant nucleons, (Np), is the same for central S+S and Pb+Pb 
collisions 158 A-GeV and equal to (h~) / (Np) = 1.93+0.14. We assume therefore the same 
value of the ratio for inelastic Pb+Pb collisions. Using the mean number of participant 
nucleons for the latter collisions calculated within the Fritiof model p2| ((Np) = 102) we 
get (h~) = 197 + 14. This leads to (J/if>)/(hr) = (1.05 + 0.11) • 10~ 6 for inelastic Pb+Pb 
collisions at 158 A- GeV. 

A similar procedure is used to calculate the (J/if))/{h~) ratio for O+Cu, O+U and 
S+U interactions. There are however two difference. The J /if) cross sections for oxygen 
and sulphur induced reactions are measured at 200 A-GeV Therefore for the compar- 
ison with the results at 158 A- GeV the measured values are scaled down by 25% according 
to the energy dependence of the J /if) multiplicity established for p+p interactions (see 
Section 2). Due to projectile-target asymmetry the Xp distribution of J /if) is expected to 
be not symmetric with respect to reflection at xf = 0. The correction for this effect is 
negelected. 
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The ratios obtained for nucleus-nucleus collisions and the corresponding ratio for N+N 
interactions at the 158 A-GeV are shown in Fig. 2 as a function of (Np). It is surprising 
that the ratio (J/ip)/ (h~) is similar for nucleon-nucleon and nucleus-nucleus interactions. 
One should however keep in mind that the ratios for nucleus-nucleus collisions may be 
underestimated by up to 30%. We repeat that this uncertainty is due to limited ac- 
cepetance of the J/ip measurement for nucleus-nucleus collisions. This systematic error 
can be reduced when the results on the rapidity or xp distributions are published. 

Finally we note that the ratio {hard process) / (h~) is expected to increase by a factor 
of about 3 when going from N+N to Pb+Pb interactions, where {hard process) denotes 
here a mean multiplicity of any process for which the cross section in A+A collisions 
increases as A 2 . 
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4 J/ip Production in p+A Interactions 



The inclusive cross section for J/ip production in p+A interactions is measured in the 
region xp > and it is usually parametrized as [Q: 

a J ^ = a (J/^)-A a{J ^\ 

where <7o(J/V ; ) and a(J/ip) are parameters fitted to the experimental data. A strong 
dependence of a{J/ip) on xp was recently measured by the E866 Collaboration |23|] at 



800 GeV. The xp distribution of J/tp decreases by a factor of about 10 from xp = to 
xp = 0.4 [21]. Thus the A-dependence of the integrated J/ip cross section in the region 



xp > is dominated by the dependence measured close to xp = 0. The values of a(J/ip) 
obtained from xp integrated data (xp > 0) or from the data close to xp = range from 



0.89 to 0.94. The results were obtained by various experiments [|1(], 27] in the 

collision energy range 200-800 GeV and they were compiled in [|J. The a(J/ip) values 
are shown in Fig. 3 as a function of \/s (filled circles). 

In order to compare the A-dependence of J / ip and h~ production the parameter a was 
fitted here to data |15] on the total multiplicity of negatively charged hadrons produced 
in p+A interactions at 200 GeV and 360 GeV. In the fit the multiplicity of proton- 
nucleon (p+N) interactions at the corresponding energy was included. This multiplicity 
was calculated as (h~) = ((h~) pp + (h~) pn )/2 |15| . Finally the a parameter fitted to 
the multiplicity data was added to the a parameter obtained by the fit to the inelastic 
cross section results (a = 0.72 ± 0.01) The obtained values of a(h~) (a(h~) = 0.88 
± 0.01 at 200 GeV and a(h~) = 0.90 ± 0.02 at 360 GeV) are shown in Fig. 4 (open 
circles). The values of a(h~) are similar to the values of a(J/ip). There is no evidence 
for any significant energy dependence both for a(h~) and a(J/ip). Similar values of the 
a parameter for hr and J/ijj production imply that the ratio (J/ip)(xf > 0)/(h~) is 
approximately independent of A for p+A interactions at high energy. We note that the 
difference in the a parameter of 0.02 (typical for the values shown in Fig. 3) results in 
10% change in the multiplicity ratio between p+N and p+Pb interactions. This can be 
compared to about 70% increase of the ratio (hard process) / (h~) expected when going 
from p+N to p+Pb interactions, where (hard process) denotes here a mean multiplicity 
of any process for which the cross section in p+A interactions increases as A (a = 1). 

The measurements of J/ip production in the backward hemisphere (xf < 0) are poor. 
However, the experimental data ]23| seems to indicate that a(J/ip) for xp < is similar to 
a(J/ip) for xp > (we note that this is not the case for pion production). This suggests 
that our conclusion concerning the similar A-dependence of h~ and J/ip production, 
based on the J/ip data from the forward hemisphere only, may remain unchanged when 
the J/ijj results in full phase space become available. 
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5 Conclusions 



The main result of this paper is that the ratio of mean multiplicities of J/i[) mesons and 
pions is similar for inelastic proton-proton, proton-nucleus and nucleus-nucleus inter- 
actions at 158 A-GeV. In our opinion this experimental observation justifies a question 
whether the generally accepted picture of J / if) creation based on the factorization theorem 
of the perturbative QCD and subsequent suppression of the J / if) yield by the interactions 
with the surrounding medium is valid. In this picture the observed scaling behaviour of 
the data, (J /if))/ (h~) « const(A), can be treated only as due to accidental cancelation of 
several large effects. 

It is obvious that the mechanism of J/if> production can not be understood without 
data on open charm creation. Published data on D and D production in p+A interactions 
are insufficient. The results on the A-dependence are summarized in Fig. 4, where 
a(D,D) is shown as a function of xp for interactions at 400 GeV [^9| and 800 GeV |30| . 
It is clear that these data do not allow to distinguish between a ~ 1, as usually assumed for 
charm production on the basis of perturbative QCD, and a ~ 0.9, the value obtained for 
pion and J/ip production. Data on open charm production in nucleus-nucleus collisions 
do not exist. It is therefore crucial for our understanding of the mechanism of charm 
creation and J / if) production to measure open charm yields in nucleus-nucleus collisions. 



Acknowledgements 

I thank Reinhard Stock for triggering publication of this work. I also thank him 
Andrzej Bialas, Leonid Frankfurt, Mark I. Gorenstein, Jorg Hiifner, Helmut Satz, Peter 
Seyboth and Herbert Strobele for discussions and comments. Claudie Gerschel and Carlos 
Lourenco helped me to understand data on J /if) production. 



7 



References 

[1] J. C. Collins, D. E. Soper and G. Sterman, in Perturbative Quantum Chromodynam- 
ics, ed. A. H. Mueller, World Scientific, Singapore, 1989. 

[2] Hard Processes in Hadronic Interactions, eds. H. Satz and X.-N. Wang, Special Issue 
of Int. J. Mod. Phys. 10 (1995) 2881. 

[3] C. Gerschel and J. Hiifner, Z. Phys. C56 (1992) 171. 



[4] C. Gerschel and J. Hiifner, Charmonium suppression in heavy-ion collisions, |hep- 
| ph/9802245] (1998). 

[5] H. Satz, Colour Deconfinement and J/ip Suppresion in High Energy Nuclear Colli- 
sions, |hep-ph/971 1289| (1997). 

[6] T. Alber et al. (NA35 Collab.), Eur. Phys. J. C2 (1998) 643, 

S. V. Afanasjev et al. (NA49 Collab.), Nucl. Phys. A610 (1996) 188c, 
H. Appelshauser et al. (NA49 Collab.), Baryon Stopping and Charged Particle Dis- 
tributions in Central Pb+Pb Collisions at 158 GeV per Nucleon, |nucl-ex/98 100141 
(1998), submitted to Phys. Rev. Lett., 
M. Gazdzicki, J. Phys. G23 (1997) 1881c. 

[7] M. Gazdzicki and M. I. Gorenstein, On the Early Stage of Nucleus-Nucleus Colli- 
sions, |hep-ph/9803462| (1998). 

[8] A. Bamberger et al, Nucl. Phys. B134 (1978) 1. 

[9] M. J. Corder et al, Phys. Lett. 98B (1981) 220. 

[10] J. Badier et al., Z. Phys. C20 (1983) 101. 

[11] E. Nagy et al, Phys. Lett. 60B (1975) 96. 

[12] C. Morel et al, Phys. Lett. B252 (1990) 505. 

[13] C. Caso et al. (Particle Data Group), Eur. Phys. J. C3 (1998) 1. 

[14] G. A. Schuler et al., Cern Report, CERN-TH. 71 70/94 (1994). 

[15] M. Gazdzicki and D. Rdhrich, Z. Phys. C65 (1995) 215. 

[16] C. Baglin et al. (NA38 Collab.), Phys. Lett. B270 (1991) 105. 

[17] M. C. Abreu et al. (NA50 Collab.), Phys. Lett. B410 (1997) 327. 

[18] A. Romano et al. (NA50 Collab.), Proceedings of XXXIII Moriond Conference, March 
22-29, 1998. 

[19] J. C. Collins and D. E. Soper, Phys. Rev. D16 (1977) 2219. 



8 



[20] V. Friese et al. (NA49 Collab.), J. Phys. G23 (1997) 1837c. 

[21] M. Anikina et al. Yad. Fiz. 38 (1983) 149. 

[22] H. Pi, Computer Physics Commun. 71 (1992) 173. 

[23] R. S. Towell et al. (E866 Collab.), Proceedings of XXXIII Moriond Conference, March 
22-29, 1998. 

[24] M. H. Schub et al, Phys. Rev. 52 (1995) 1307. 

[25] D. M. Aide et al. (E772 Collab.), Phys. Rev. Lett. 66 (1991) 133. 

[26] M. I. Leitch et al. (E789 Collab.), Phys. Rev. D52 (1995) 4251. 

[27] F. Fleure (NA38 and NA51 Collab.), Ph. D. Thesis, Ecole Polytechnique, Palaiseau 
(1997). 

[28] A. S. Carroll et al, Phys. Lett. 80B (1979) 319. 
[29] M. E. Duffy et al. (E613 Collab.), Phys. Rev. Lett. 55 (1984) 1816. 
[30] M. J. Leitch et al. (E789 Collab.), Phys. Rev. Lett. 72 (1994) 2542. 
[31] E. Cobbaert et al. (WA78 Collab.), Phys. Lett. B191 (1987) 456. 
[32] C. A. Alves et al. (E769 Collab.), Phys. Rev. Lett. 69 (1992) 3147. 
[33] M. Adamovich et al. (WA82 Collab.), Phys. Lett. B284 (1992) 453. 
[34] M. Adamovich et al. (Beatrice Collab.), Nucl. Phys. B495 (1997) 3. 



9 



Table 1 The results on mean multiplicity of J/ip mesons produced in p+p interactions. 



[GeV] 


(J/tp) ■ 10 6 


Reference 


6.8 

8.7 

19.4 

24.3 

52.0 


0.021 ± 0.006 
0.075 ± 0.037 
3.8 ± 0.3 
4.6 ± 0.8 
19.7 ± 8.7 
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Figure 1: The multiplicity of J/ip mesons produced in p+p interactions as a function 
of the collision energy. The filled circles indicated measured data. The open circle shows 
the estimated multiplicity at y/s = 17.3 GeV. 
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Figure 2: The ratio of the mean multiplicities of J/ip mesons and negatively charged 
hadrons for inelastic nucleon-nucleon (square) and inelastic O+Cu, O+U, S+U and 
Pb+Pb (circles) interactions at 158 A-GeV plotted as a function of the mean number 
of participant nucleons. For clarity the N+N point is shifted from (N P ) = 2 to (N P ) = 5. 
The dashed line indicates the mean value of the ratio. 
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Figure 3: Comparison between a(J/ip) (filled circles) and a(h~) (open circles) for p+A 
interactions in the energy range 200-800 GeV. The dotted line shows the value a = 1 
characteristic for the A-dependence of total charm cross section obtained in models based 
on the perturbative QCD. The dashed line indicates the value a = 0.9 measured for pion 
production in full phase space. 
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Figure 4: 
and 7T+A interactions 



Dependence of a(D,D) on for p+A interactions [29, |3(J (filled symbols) 
,§1|| (open symbols) at 250-800 GeV. Circles indicate 
results obtained by reconstruction of D and D decays. Squares indicate data obtained 
by the analysis of prompt single leptons or neutrinos. The results for which the (xp) is 
not given are plotted at the lower edge of the acceptance region. The dotted line shows 
the value a = 1 characteristic for the A-dependence of total charm cross section obtained 
in models based on the perturbative QCD. The dashed line indicates the value a = 0.9 
measured for pion production in full phase space. 
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